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Abstract 

We present an extension of the density-functional theory (DFT) 
formalism for lattice gases to systems with internal degrees of free- 
dom. In order to test approximations commonly used in DFT ap- 
proaches, we investigate the statics and dynamics of occupation 
(density) profiles in the one-dimensional Potts model. In particu- 
lar, by taking the exact functional for this model we can directly 
evaluate the quality of the local equilibrium approximation used 
in time-dependent density- functional theory (TDFT). Excellent 
agreement is found in comparison with Monte Carlo simulations. 
Finally, principal limitations of TDFT are demonstrated. 
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1 Introduction 



Density functional theory (DFT) and its time-dependent variant (TDFT) are powerful 
methods to derive phase diagrams and the kinetics of phase transformations |H] in 
condensed matter systems, in particular in the presence of confinement effects. Particularly 
useful are these theories for lattice systems (UEj, allowing us to deal with the discrete nature 
of structures encountered e.g. in the description of metallic alloys, adsorbate layers, and 
complex pattern formation on atomic scales. 

In the case where many ordered phases can coexist one is lead to include internal degrees 
of freedom into the DFT approach. As a prominent model we will consider in this work the 
(/-state Potts model, which has a (/-fold degenerate ground state The non-equilibrium 
dynamics of that model reflects the coarsening of domains following a quench from the 
disordered homogeneous phase to a system with long-range order seen in binary alloys, 
liquid crystals, magnetic bubbles, Langmuir films and soap bubbles [7j. The Potts model 
(in the limit q — ► 1) is isomorphic to a site-bond percolation problem [8] EJ and for q = 2 
it corresponds to the Ising model jlUj . There are interesting experimental realizations for 
q = 3 (e.g. Kr on Graphite 9 = 4 and q = oo |12j (froths and metallic grains). 

In this paper, based on the exact density functional for the one-dimensional (/-state 
Potts model [lEl, we investigate the quality of various approximations often employed in 
DFT and TDFT. After presenting a general scheme for treating systems with internal 
degrees of freedom, we first consider the mean spherical approximation (MSA) for the 
equilibrium properties, and show how it compares with a mean-field approach and the 
exact solution. We then derive the TDFT ^1] for the one-dimensional Potts model and, 
as an application, study the smoothening of an initial sharp-kink density profile. Let 
us note that this TDFT differs in many respects from the Runge-Gross theory ^S] for 
electronic systems in a time-varying external field, one major difference being that the 
present approach assumes stochastic, overdamped dynamics. The basic techniques to set 
up a TDFT for classical systems are reviewed in |16( Oj . By comparison with Monte-Carlo 
simulations we show that the TDFT provides a significant improvement over kinetic mean 
field theory. Moreover, since in our case the TDFT is based on an exact density functional, 
the differences between the simulated and DFT results allow us to perform a specific test 
of the local equilibrium approximation used in the TDFT. Finally, we discuss principal 
limitations of the TDFT, where only the density profile is used to specify a non-equilibrium 
state of the system. This amounts to an incorrect account of correlations. 

2 DFT for Lattice Gases with Internal Degrees of Freedom 

We consider a lattice gas where each site i is either vacant {xf = for all a) or singly 
occupied by a particle in one of q internal states (note that this corresponds to a generalized 
Potts model with q + 1 states). If state a € {1, . . . , q} is realized, xf = 1 and = for 
j3 ^ a. Occupation numbers therefore satisfy xfxf = xf5 aj 3. The Hamiltonian including 
two particle interactions and site energies ef due to an external potential is given by 



The formal steps of DFT for fluids or lattice gases without internal degrees of freedom can 
be carried over to the case considered here. One arrives at a variational principle based on 
the functional for the grand canonical potential 
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— Ha, fJ-a being the chemical potential fixing the mean total occupation p a 
of state a, p = {pf } and pf = (xf) are the average occupation numbers. The free energy 
functional is decomposed into an 'ideal' part describing a non-interacting lattice gas, 
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and an excess part -F e x[p] due to interactions (for convenience we set k^T = 1). The 
equilibrium occupation is then obtained by minimizing J)[p] with respect to the pf. The 
corresponding equations dQ/dpf = determining the equilibrium profile are called struc- 
ture equations. The minimum value of f2 [p] is the grand-canonical potential at equilibrium. 

Higher order derivatives of £1 [p] with respect to the pf taken at the equilibrium profile 
yield a hierarchy of direct correlation functions |2j. In particular, by differentiating the 
structure equation, the inhomogeneous Ornstein-Zernike equation 
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can be derived [I], which relates the direct correlation function c 
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degrees of freedom is recovered. 
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1 the known result for systems without internal 



3 Application to the ID Potts model 

In the following, standard approximation schemes used in DFT are tested based on the 
one-dimensional Hamiltonian (1 < i < M) 

# =E<^m + E e >"> ( 5 ) 
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which for vf 13 = V5 a p reduces to the standard Potts model j^j. 



3.1 Exact density functional 

For one-dimensional lattice gases with short range interactions, a general scheme for deriv- 
ing exact density functionals based on Markov chains has recently been developed in |13j . 
This approach in particular provides an exact functional for the Potts model, which has 
been derived earlier by Percus ^7j. The functional can be written as 
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where the correlators 

rf = (xjLxxf > (7) 

have to be expressed by the mean occupation numbers {pi}- This is achieved by solving 
the correlator equations 
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A simpler expression for the functional is obtained if the vacancies are considered as an 
additional Potts state with index a = (and vf?° = = 0), whereby the lengthy entropic 

part reduces to YliLi Yla 8=o ^°s(^i/pf)- in the fully occupied case i^2, a xf = 1) the 
solution of eqs. (JSJ) requires a careful limiting procedure by letting the vacancy concentration 
go to zero. 

The structure equations read 
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By solving this equation numerically we obtain the exact density profiles. (For the example 
of eq. (15), results are shown in Fig. 1.) 

3.2 Mean spherical approximation (MSA) 

In the following the excess density functional will be expanded around a homogeneous 
reference state, pf = p a |1H] . For simplicity, we here consider the standard Potts case 
= VS a p, and set p a = p for all a. 

Defining Apf = pf — p, AO[p] = f2[p] — f2[{p}], etc., we can write 

An[p] = AF id [p] + AF cx [p] + £ If Apf . (10) 

When retaining only terms up to second order in the excess free energy functional, we find 
AF cx [p] = -Y J c {1 Hp)Apf - \ c^(\i-j\,p)ApfApl (11) 

i,a i,a,j,(3 



where c^(p) can be subsumed into the chemical potential. For our choice of interactions 
with translational symmetry, g?? = g a p{\i — j\), and the direct correlation function ^ 



can be split into two parts, 
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where in the MSA approximation, 
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and g a p(l = 0) = 
Zernike equation 



0. The two unknowns ci(0) and 02(0) are found using the Ornstein- 
P, which yields 
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As an example, we now consider the special case of zero external energies for the Potts 
states a = 2, . . . ,q and alternating external energies for the Potts state a = 1, 



e, % even and a 
0, otherwise 
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and assume fi a = fx independent of a. Due to the symmetry of this external potential, 
each occupation number pf is equal to one of the four representatives Pq,Pq,p\,Pi, and the 
mean particle density per site is 
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We consider the functional 
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with /i e ff = jU— e/2 + c^ 1 ^ (p). The structure equations are obtained by setting the derivatives 
of \£[p] with respect to p\, pfi, p\, p\ equal to zero. These equations are solved numerically 
subject to a fixed n in eq. lfT7)|) (which is equivalent to d^/dfi e ff = 0). Ordinary mean field 
theory can be recovered from the functional by setting ci(0) = 02(0) = 0. 

In Fig. 1 the occupation numbers for the MSA approximation and the mean field theory 
are compared with the exact results based on eq. @. Without interactions (V = 0) the 
external potential yields mean occupation numbers p$ < p\ = pf < Po, p$ = e~ e Po- When 
a repulsive interaction V > is switched on, the fact that p\ > p\ induces an increasing 
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Figure 1: Comparison of the occupation numbers Pq, p\ and p\ in dependence of interaction 
strength V for two q values in the MSA approximation (figures on left side, dashed lines) 
and mean field approximation (figures on right side, dashed lines) with the exact solution 
(solid lines) for n = 0.8 and e = 1. The corresponding MSA and mean field results for p$ 
(not shown) are always close to the exact solution. 



occupation difference. For larger q this effect becomes less pronounced due to the increasing 
contribution of the entropy to the free energy per site. 

As expected, the MSA is an improvement over the simple mean field approximation 
especially for higher values of the interaction parameter and high overall densities n. More- 
over, the quality of the MF-results can be shown to improve when e is decreased. A notable 
result is the improvement of the two approximations for larger q, reflecting the fact that 
mean-field descriptions should become exact in the limit q — > oo. 



3.3 Kinetics of density profiles 



So far we have considered density profiles at equilibrium. In order to account for the time 
evolution of non-equilibrium profiles we will make use of the local equilibrium approximation 
for the probability distribution P(x, t) to find a configuration x = {xf} at time t. In this 
time-dependent density functional theory (TDFT) [HI El, the deviations of P(x, t) from the 
Boltzmann equilibrium distribution are described by a one-particle time-dependent effective 
potential hf(t), 
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The effective potential is the unique potential, which yields the instantaneous density profile 
pf(t) according to the equilibrium DFT. Accordingly, the unknown field h in eq. (|18|) can 
be determined by the "structure equation" hf(t) = — dVt\p\/dpf with f2[p] from eq. (J2J- It 
is then clear that in this approach all equilibrium relations between occupational correlators 
(xfxj . . .) and the density p also apply at each time instant to the non-equilibrium situation 
[in particular eq. ©]• 

To test the quality of the TDFT, we consider a lattice with all sites being occupied 
(n = 1) and a nonconserved dynamics, where a given Potts state a on lattice site i can 
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change to any other state (3 ^ a with a rate tf"^(x) that depends on the current state x 
due to interactions. From the master-equation describing this stochastic process, we derive 
the equation of motion for the occupation profile, 
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where (■ ■ -)t denotes an average with respect to the non-equilibrium distribution P(x, t). 
For the rates we choose a generalized Glauber form, 



w 



a/3 



(x) 



1 — tanh 



- xtxxf )V/2 



1 — tanh 



(xf +1 xf-xf +1 xf)V/2\) , 



(20) 



which satisfies the condition of detailed balance |20| . 

When substituting eq. (j20j) in ()19jl further evaluation is made possible by noting that the 
factor in front of V/2 in the argument of tanh is always —1,0 or 1, so that it can be taken 
out from the tanh function. With the help of the Markov property, three-point correlators 
can exactly be reduced to two-point correlators, e.g. (xf_ x x^ xj +l ) = T^J-^/p^. After 
some algebra we then find 
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Equations (|21jl together with the correlator equations (jSJ) form a complete set of equa- 
tions for the time evolution of density profiles and can be solved numerically for a given 
initial condition. 

As an example, we consider an initially sharp kink profile at time t = 0, where the left 
part of the system is in Potts state a = 1 and the right part in Potts state a = q with fixed 
boundary sites i = 1 and M in Potts states a = 1 and q, i.e. x\ = S a> x and = S a>q . In 
Figure we compare the time evolution of the profile with the results from continuous- 
time Monte-Carlo simulations. The agreement is excellent for all times until for large 
times the TDFT solution and the simulations both yield the correct equilibrium profile. 
In order to obtain this agreement, it was necessary to adjust the time scale according 
to tTDFT = 0.85 tyic- The significant improvement over a simple mean-field treatment 
corresponding to factorization of all correlators in ea.(|21j). i.e. T? 13 = p?_iP?, etc., can be 
seen by comparing Fig. [2ji with Fig. [2j3. The mean-field approximation is insufficient from 
the beginning and does not provide the correct equilibrium profile for long times. 



3.4 Basic Limitations of TDFT 

Since in the kinetic equation of the TDFT only the densities (mean occupation numbers) 
enter, the TDFT can not distinguish between states with the same density profile but 
different correlations. As a consequence, the time evolution of a non-equilibrium state 
with the equilibrium density profile cannot be captured. This issue may become important 
when considering memory effects in systems with slow relaxations, as for example glassy 
systems. The observables used to characterize the thermodynamic states of such systems 
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Figure 2: (a) Comparison of the evolution of an initial kink profile calculated using TDFT 
(lines) with Monte Carlo simulations (points) for a system with 32 particles, q = 4 and 
V = —2. Boundary conditions are pf = 6 a> i and pfy = 5 a ^. Results are shown for 
t = 1,2,4,8 and the equilibrium density (thick lines). Results for a = 4 (not shown) are 
mirror images of those for a = 1 with respect to the center site i = 16. (b) Comparison of 
mean field theory (dotted lines) with the simulations (solid lines) for the same system as 
shown in (a). Equilibrium solutions are indicated with thick lines. 

(as e.g. density) may have equal values initially, but very different time evolutions depending 
on the systems' history. 

To illustrate this principal failure of the TDFT, we generate states of the Potts model 
with the same mean occupation numbers but different correlations for the same system 
as in the previous chapter. To achieve this, we first determine the equilibrium profile 
pf for a system with interactions from the structure equations ©• This profile is also 
the equilibrium profile for a system without interactions (V = 0) and external potential 
ef = — logp^- Hence by taking the equilibrium state of the non- interacting system as initial 
non-equilibrium state for the interacting system, we can study by Monte-Carlo simulations 
the time development of a density profile that according to TDFT cannot change. 

Figure |3] shows the effect for V = — 4 and q = 4. Although the mean occupation 
numbers at time t = are the same as in the equilibrium state reached for t ~ 1024, there 
is a pronounced change in the occupation profile at intermediate times. These changes are 
due to the fact that in the initial state the correlations between occupation numbers are not 
the equilibrium ones. The configurations with large statistical weights in the initial state 
exhibit stronger short-wavelength fluctuations and in order for the correlations to build up, 
the system has to pass through intermediate states with a non-equilibrium profile. 
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Figure 3: Evolution of a density profile starting from the equilibrium density, but with 
non-equilibrium correlators (see text) for V = —4 and q = 4. The system quickly reaches 
the state with maximum deviation from the equilibrium density at t ~ 1 and then relaxes 
to the original equilibrium density profile. Only the left half of the system is shown. 
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